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Abstract 

We use field theoretic renormalization group methods to study the 
critical behavior of a recently proposed Langevin equation for driven 
lattice gases under infinitely fast drive. We perform an expansion 
around the upper critical dimension, d c = 4, and obtain the critical 
exponents to one loop order. The main features of the two loop calcu- 
lation are then outlined. The renormalized theory is shown to exhibit 
a behavior different from the standard field theory for the DLG, i.e. 
it is not mean field like. 

Since it was first introduced by Katz et al. 0, the driven lattice gas model 
(DLG hereafter) has attracted considerable interest M [jj. Being one of the 
simplest archetypes of non-equilibrium model its study may contribute to 
pave the way for an understanding of out of equilibrium systems. The DLG 
consists of a periodic regular lattice on which nearest-neighbor particle-hole 
exchanges are performed. The hopping rate is determined by the energetics 
of the Ising Hamiltonian H, the coupling to a thermal bath at temperature 
T, and an external uniform driving field E pointing along a specific lattice 
axis. In particular, the hoping rate depends on [(AH + £E)/T], where AH 
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is the energy variation which would be caused by the tried configuration 
change, E = |E|, and I = 1 (-1) for jumps along (against) E and other- 
wise (see [0, |3| for a detailed description). The DLG exhibits a continuous 
phase transition from a disordered state at high T to a stripe like ordered 
state at sufficiently low T 0, ||. The nature and properties of this transition 
have been largely studied in recent years. A new general Langevin equation 
has been proposed, capturing the physics of the DLG at the critical point 
||]. This Langevin equation predicts different critical behavior for the cases 
< E < oo and E = oo (in which particles cannot jump against the field) 
respectively [B, [5]. That is, the point E = oo behaves as a sort of tricritical 
point in the parameter space. For finite values of the driving field E the 
Langevin equation previously proposed by Janssen and Schmittmann is re- 
covered H . It is the purpose of this paper to investigate the critical behavior 
of the DLG for E = oo in order to determine explicitly the differences with 
the < E < oo case. The new Langevin equation reads M, Q 



A||Aj> - Ai0 + rA ± + |Aj 



+V^V ± -£ ± + y| V,|£||, (1) 

where Vn (VjJ is the gradient operator in the direction parallel (perpendic- 
ular) to the electric field, and the noise satisfies 

(£CM)> = o, 

(V.€(x,f)-V / .*(*',0> = -V 2 5(x-x')<5(t-t')- (2) 

A very similar equation has been proposed to describe the Freedericksz tran- 
sition in nematic liquids, and general asymmetric two-dimensional pattern 
formation 0. 

In order to renormalize this equation, following standard methods || , let 
us introduce a Martin-Siggia-Rose response field and recast Eq. ([!]) as a 
dynamical functional M, the associated action of which is 



£U 6) = / d d xdt 



e o, 



d t - y(-A||A_L - Ai + rA ± )j0 - ^0A ± 
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Figure 1: Elements of perturbation theory: the response and correlation 
propagators and the four-point vertex. 0-legs are indicated by a wiggle. 



The free propagators are 
G° 2 (k,u;) 



u 2 + {^) 2 k\{k 2 + r) 2 ' 
<ZM = i^W(^7Y (4) 

and the vertex is: —eog/12 k\. These elements can be represented diagram- 
matically as in Figure 1 (wavy legs symbolize response fields; straight lines 
stand for density fields). 

In order to renormalize the theory, one has to look for the primitive diver- 
gences in a perturbation expansion. If r An denotes a one-particle irreducible 
vertex function with h external 0-legs and n external 0-legs, only Tn and 
Ti3 are found to possess primitive divergences. The Feynman diagrams con- 
tributing to these vertex functions are topologically identical to model B 
graphs M. However, the bare correlation and response propagators that fol- 
low from Eq. (|3]) are anisotropic, in contrast to their counterparts in model 
B§. 

To one loop in e = 4 — d, the ultraviolet divergences in r n and T 13 lead to 
the renormalization of r and u, the latter being the dimensionless coupling 
constant u = A e r~ £ l 2 g. A £ is a numerical factor to be defined below. We 
define renormalized parameters r« and Ur by = Z T r and Ur ; = Z u u. 
Given that the leftmost diagram in Figure 2 does not depend on external 
moments or frequencies, the derivatives of r n with respect to them vanish, 
and no extra (field) renormalizations are therefore required. The Z factors 
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Figure 2: One loop diagrams contributing to r n (left) and T 13 (right), 
are determined by the following normalization conditions 

a r iJ| e o 

"k*_ L n\NP = ~^ T R' 

d ki r«\ NP = ^r^W (5) 

A convenient choice for the normalization point NP is hi = Ui = and 
r = /i 2 , where fi is an arbitrary momentum scale. To one loop, we find 

T n (w,k) = iuj + ^kKk 2 + r) + D u 

T 13 (w,k) = ^k 2 ± g + D 2 , (6) 

where Di(D 2 ) corresponds to the algebraic expression of the left(right) dia- 
gram in Figure |2|. A calculation in dimensional regularisation |10|, || yields 



D 1 
D 2 



5^e 2 rl 6/2 



1287T 2 



% 2 e , 2 T £/2 



1287T 2 



(7) 



After setting A £ = 5/647T 2 , one obtains 



9^ ^iz\np 



e 5ge r 1 e/2 
— r 



e 



1287T 2 



1 + 



5t s/2 
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Mn 2 e_ 



e 



1 + 



u 



2 £ 



e. 



(8) 
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Figure 3: Two loop contribution to r n . 



which entails 



u „, 2 



Z T = 1 + - + 0(V 

£ 

Z u = l-- + 0(u 2 ). (9) 

The renormalization group equation obtained after requiring invariance of the 
bare irreducible vertex functions upon changes on the normalization point 
reads 

[fid, + pd UR + cd TR \rt = o, (io) 

where the renormalization group functions are defined in the usual way: 
f3{u R ) = fid^UR, and ((u R ) = /i<9 M (ln r^). A straightforward calculation 
then leads to 

(3(u R ) = -eu R + u 2 R + 0(u R ), 

CM = 2-u R + 0(u 2 R ), (11) 

from which one can determine the location and stability of the fixed points. 
To this order, apart from the trivial mean-field result u* R = 0, a nontrivial 
fixed point u* R = e emerges. This fixed point controls the critical behavior of 
the theory below four dimensions. 

Now we proceed with the calculation of the associated critical exponents. 
We first note that, as indicated above, no renormalization of the fields 0, 
4> has been required. Therefore, in particular, the anomalous dimension of 
vanishes, i.e., r\ = 0. Concerning the exponent z/, which controls the 
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divergence of the correlation length with temperature, we simply have v = 
((u^)' 1 = 1/2 + e/4 + 0(e 2 ). This is to be compared with v = 1/2, the 
value obtained by Janssen and Schmittmann for E = oo ||. This result 
demonstrates that the continuous version Eq. (|l|)of the DLG with E = oo is 
not mean-field like but characterizes a universality class other than the one 
in ||. Since there are no dangerous irrelevant operators in Eq. ([|) standard 
scaling laws apply (contrary to the case in |J). Therefore the exponents 
are related to each other and estimating rj and v is sufficient to deduce all 
the other exponents. For instance, the order parameter exponent (3 can be 
written as (3 = |(d - 2 + rj) |TtJ, and we have = 1/2 + 0{e 2 ). 

The previous results concern the one loop approximation. The two loop 
calculation presents an important new feature, namely, that the scaling be- 
comes anisotropic. The detailed two loop calculation involves the evaluation 
of many non trivial integrals and will be presented elsewhere. We do not 
present it here, but let us stress that a simple glance to some of the contribut- 
ing integrals permits us to extract far reaching consequences. In fact, Figure 
3 reveals that, contrary to what happens to one loop order, the graph con- 
tributing to Tn to two loops depends on external frequencies and momenta. 
Therefore, eo and the fields need to be renormalized to heal these new diver- 
gences. Moreover, this diagram depends on the parallel and perpendicular 
components of the external momenta in an asymmetric way; consequently, 
anisotropic exponents emerge. On the other hand, after including two loop 
diagrams, corrections to the (3 = 1/2 mean field value may appear. Finally, 
we remark that to two loops usual scaling laws may no longer be valid; in 
particular, the exponent v will split up in u± and v\\. 

Summing up, we have performed the renormalization of the field theory 
in [|| for the DLG under an infinitely large driving field. The renormal- 
ization procedure yields results essentially different from those for a finite 
field. In particular, corrections to mean field are observed explicitly in the 
one loop approximation for the exponent v. Anisotropic exponents and a 
non-mean-field exponent (3 are predicted from simple arguments based on 
the analysis up to two loop diagrams. This calls for extensive computational 
simulations to observe numerically the difference between finite and infinite 
driving cases. Acknowledgements It is a pleasure to acknowledge J. Marro, 
J. L. Lebowitz and P.L. Garrido for useful discussions, and E. Hernandez- 
Garci'a for pointing out reference @ to us. This work has been partially 
supported by the European Network Contract ERBFMRXCT980183 and by 
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